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Abstract The coherence length of the thermal electromagnetic field near a planar 
surface has a minimum value related to the nonlocal dielectric response of the 
material. We perform two model calculations of the electric energy density and the 
field's degree of spatial coherence. Above a polar crystal, the lattice constant gives 
the minimum coherence length. It also gives the upper limit to the near field energy 
density, cutting off its 1/z 3 divergence. Near an electron plasma described by the 
semiclassical Lindhard dielectric function, the corresponding length scale is fixed 
by plasma screening to the Thomas-Fermi length. The electron mean free path, 
however, sets a larger scale where significant deviations from the local description 
are visible. 

PACS: 42. 25. Kb Coherence - 07.79.Fc Near-field scanning optical microscopes 
- 44.40.+a Thermal radiation - 78.20.-e Optical properties of bulk materials and 
thin films 



1 Introduction 

Thermal electromagnetic radiation in vacuum, as described by the celebrated black- 
dody spectrum discovered by Max Planck 1 1 1, is usually taken as a typical example 
of incoherent radiation. This is not quite true, however: if the radiation is detected 
at a given frequency, it is spatially coherent on a scale set by the wavelength 82131 . 
When one approaches a macroscopic object, the radiation spectrum and its coher- 
ence is profoundly changed, depending on the properties of the object. For exam- 
ple, if the object supports resonant modes like surface plasmon polaritons, the field 
is coherent across the propagation length of these modes |4|. The opposite case is 
possible as well: the coherence length becomes comparable to the observation dis- 
tance, much smaller than the wavelength, close to an absorbing object with a local 
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dielectric function [5 1. It has been suggested already by Rytov and colleagues that 
this behaviour is an artefact because at some small scale, nonlocal effects must 
come into play |2|. This is what we discuss in this paper in a semi-quantitative 
way. We use two models for nonlocal dielectric functions and identify the scale 
for the field's coherence length using explicit asymptotic formulas. A nonlocal di- 
electric response is of primary importance for semiconductor quantum wells, see 
for example Ref.|6|, but the issue of spatial coherence has not been analyzed in 
this context, to our knowledge. 

We focus on the spatial coherence of the electromagnetic field at nanometer 
distance in the vacuum (medium 1) above a solid material (medium 2). We chose a 
planar geometry which is sufficiently simple to allow for an essentially analytical 
description, thus avoiding the need for extensive numerics. On the other hand, 
many effects have been discussed in this setting: the fluorescence rate of molecules 
near metals and thin films |7 1, scanning near- field microscopy of sub-wavelength 
objects deposited on a substrate |8|, the momentum exchange between a tip and 
a sample (related to the Casimir force, see, e.g., |9|) and the energy exchange 
between a tip and a sample 11011 111 211 31 . 

2 Basic notation 

2.1 Field correlations 

The spatial coherence of the electric field is determined by the two-point expecta- 
tion value 1141 



where the average is taken in a stationary statistical ensemble (thermal equilibrium 
in the simplest case). We focus in the following on the cross-correlation spectrum 
£ij (ri , r 2 ; oS) and a frequency in the infrared to visible range. Far from any sources 
and in global equilibrium, the corresponding wavelength A = 2nc/u) sets the scale 
for the field's spatial coherence length: the cross-correlations tend to zero if the 
distance |ri — r2 1 exceeds A. In the vicinity of a source, the coherence length ^ co h 
significantly differs from A, as Henkel and co-workers have shown previously |5 1, 
and it changes with the observation point. 

The spectrally resolved electric energy density is given by the trace 



and its value in thermal equilibrium allows to define an electric, local density of 
states, as discussed in more detail by Joulain and co-workers 1151 . The normalized 
tensor 
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to be considered below, allows to introduce a spatial degree of coherence. In the 
following, we call a "coherence function" both, £jj(ri, T2\ijS) and Eq.(|3jl- Defini- 
tions for a degree of polarization based on this 3x3 matrix (with ri = r-2) have 
been put forward as well, see 1161171 . For the sake of simplicity, we suppress the 
frequency arguments in the following. 

2.2 Planar surface with local response 

In a previous paper, Henkel and co-workers have shown that in the vacuum above 
a planar dielectric surface at temperature T, described by a local permittivity £2, 
the spatial coherence function is of the form |5 1 (see also II 81191 ) 



where 0(uj,T) — hu/(e hLO / kT — 1). We assume that the field is observed in 
vacuum (relative permittivity e\ = 1). The surface is given by z — 0. We have 
chosen the x-axis such that ri — r"2 lies in the xz-plane and p = x\ — X2- The 
quantity f 2 =p 2 + (z\ + Z2) 2 is the (squared) distance between ri and the image 
point of r"2 across the interface. 

Eq.0 applies to leading order when both distances z\, z^ are much smaller 
than the wavelength A; for other regimes and higher order corrections, see Ref.|5 1 
and, at p = 0, Ref.[20|. In the following, we focus on the correlation function at a 
constant height z = z\ = Z2 and discuss its dependence on the lateral separation 
p; note that p can be positive or negative. The normalized coherence function @ 
is qualitatively similar to a Lorentzian: the yy-component, for example, follows 
a law ~ [4z 2 + p 2 ]" 3 / 2 . The spatial coherence length is thus equal to z, and 
decreases without apparent limitation as the surface is approached. The electric 
energy density derived from @ diverges like 1/z 3 : 



Both points have been noted by Rytov and co-workers Q, who have also argued 
that this unphysical result is due to the assumption of a local dielectric response 
down to the smallest scales. A cutoff would occur naturally in a non-local treat- 
ment or taking into account the atomistic structure of the material. This is what 
we show here in detail, using two different model calculations. Doing this, we also 
provide a basis for the phenomenological cutoff introduced recently by Kittel and 
co-workers 1 13 1 in the context of heat transfer from a hot, sharp tip into a cold, 
planar substrate. 

2.3 Overview 

We will use two models to calculate the coherence function. In both, we focus, 
as mentioned before, on the fields near a planar surface and compute the field 
correlations in the vacuum above it, at sub-wavelength distances. 




(4) 



u E {z) = [0(^T)/(87rz 3 )]Im[(£ 2 - l)/(e 2 + 1)]. 



(5) 
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The first model is based on the fluctuation electrodynamics introduced by Ry- 
tov and co-workers |2j where the sources of the field are described by fluctuating 
polarization currents below the surface. This approach relies on the fluctuation- 
dissipation theorem that links the spectrum of the polarization current to the di- 
electric function of the material below the surface. We argue that the currents are 
spatially correlated on a scale equal to or larger than the material's microscopic lat- 
tice constant. We then show that the radiation generated outside the surface shows 
a minimum coherence length given by this scale; this cuts off the divergences ap- 
pearing in a local description of the material, as noted in Refs.|2 13 j. This model 
can be applied to polar ionic crystals in the frequency domain where the dielec- 
tric response is dominated by phonon-polariton resonances. It can also cover a 
non-equilibrium situation where the surface is heated to a different temperature or 
shows weak temperature gradients I21 I22I . 

The second model describes the dielectric response of an electron plasma and 
applies to the plasmon-polariton resonances occurring in metals. We use here di- 
rectly the fluctuation-dissipation theorem for the electric field [23 24 1, restricting 
ourselves to a field-matter system in "global equilibrium". The coherence func- 
tion is determined by reflection coefficients from the surface for which we take the 
Lindhard form, taking into account the non-local response of the electron plasma. 
It is shown that the field's coherence length is limited by the Thomas-Fermi screen- 
ing length, but significant deviations from the local description occur already on 
the (typically larger) scale of the electron mean free path. 



3 Polar crystal 

3.1 Current correlations 

We assume here that the fluctuating currents that generate the radiation field, are 
correlated below a certain distance /. Above this distance, the medium response 
can be considered as local. A lower limit for I is certainly the lattice period a: 
at scales smaller than a, the concept of a continuous material characterized by a 
dielectric constant does not make sense any more. 

In this situation, the cross correlation spectrum of the fluctuating currents, 
as given by the fluctuation-dissipation theorem, is no longer delta-correlated in 
space. We choose here to smoothen the spatial delta fonction into a gaussian. The 
fluctuation-dissipation theorem for the currents thus takes the form 

- (ri -r 2 ) 2 /; 2 

01(n,^(r 2 ,^)) = 2^ £o Im[ £ (r)] &(ui,T)5kiS(ui - a/), (6) 

where r = |(ri + ra). The gaussian form for the spatial smoothing is chosen 
for convenience; another functional dependence, e.g. the model put forward by 
Kliewer and Fuchs |25 1, will lead to qualitatively similar results. 



Electromagnetic field correlations near a surface with a nonlocal optical response 



5 



3.2 Transmitted field 



We then write the cross correlation spectrum for the electric field in terms of Green 
functions and the currents. We use the convention 



Ei(r,u) = ip, uj / d 3 r' G lfc (r, r'; w)j fe (r', 

« J. 



(7) 



To proceed further in the calculation, the Green function is written as a Weyl plane 
wave expansion ([26 1 and appendix). In the present case, the Green function relates 
the current on one side of an interface to the electric field on the other side of the 
interface. It depends on the Fresnel transmission coefficients through this interface. 

Using (|6j and integrating over the half-space filled with the dielectric, one 
obtains 

v /" 27r d0 [°° KReMdK 

xe -^pco see -2I m ( 7l ), e -K 2 / 2 /4 e -Re( 72 ) 2 i 2 /4 5 * fc(K)5jfc(K) 

In the preceding equation, the wavenumberin the medium i = 1, 2 is kj = (K, ji) 
where K = K cos 9 e x + K sin 6 e y and 7? = e, (w/c) 2 — k 2 . The tensor gij (K) 
is related to the Green tensor and defined in the Appendix. 

The cross-spectral correlation function depends on four characteristic lengths: 
the wavelength A, the distance to the interface z, the locality distance I and the 
separation p between the field points. The latter is the variable considered in our 
problem. At the wavelengths we work with, we always have I <C A. When z is 
larger than A (in the far field), the factor e~ 2Im ( 71 ) 2 actually limits the integration 
over K to < K < ui/c, i.e., to propagating waves. The cross-spectral correlation 
function, in this regime, drops to when p exceeds A/2, as in the blackbody ra- 
diation field. In the intermediate regime I « z « A, the integral is dominated by 
the range w/c«A'«l/I, where the exponentials containing I are close to unity. 
Hence, the results of Ref.|5| are recovered. Finally, when z -C I, e _2Im ( 71 ) 2 and 
e -Ro(7 2 ) i /4 ^Qjjj a pp roac h unity in the relevant range |^/e2|cj/c -C K < l/l. 
This is the regime we discuss in more detail in the following. 

We note in passing that we use our calculation is based on the solution to 
the transmission problem valid for a local medium. Actually, this solution applies 
when the wave vector K is smaller than 1 / / when the medium can be described 
as homogeneous. But from l|8) one sees that whatever the values of z, there is 
anyway a cut-off in the integration over K at approximately l/l. Therefore, one 
might consider that the local expression of the Fresnel coefficients remains valid. 
We believe that our model, even if it not rigorously accurate, is useful in view of 
the insight one gains from the analytic result. 
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Fig. 1 Normalized electric energy density above a surface of silicon carbide vs. the distance 
z to the surface. The electric energy density is normalized to the electric energy density in 
the far field. The locality scale is taken as / = lnm. The SiC permittivity is described by an 
oscillator model in the visible-infrared part of the spectrum 1281 . 



3.3 Asymptotics and discussion 

Using the limit of g^j (K) for large K, we obtain from (|8j the following asymptotic 
expression for the cross spectral correlation tensor 



£ 3 -(n,r 2 )« — — — — (9) 



86>(m,T)lm(e 2 ) 
e ncj\e 2 + l\ 2 l 3 

x ^[M 3/2 + |^M 5/2 ] 

V 2fe-" 2 /i 2 o 

where M 3/2 = M(|, 1, and M 5/2 = M (§, 3, -^), and M (a, 6, z) is the 
confluent hypergeometric function |27|. When p <C i, M 3 / 2 and M 5 / 2 both ap- 
proach unity. Putting p = in the cross-spectral correlation tensor and taking the 
trace, we get the electric energy density versus z: 

z«l: u E (z) = af^Im^l. (10) 

It appears (see Fig[0 that it saturates at short z to a quantity that only depends on I 
as 1/Z 3 : the non-locality scale / thus sets the ultimate length below which the field 
properties are "frozen" to their value for 
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Fig. 2 Normalized spatial coherence function vs. lateral separation p in units of the nonlo- 
cality scale I. The nonzero components in Eq.^9} are plotted and normalized to the trace of 
the coherence tensor. 

When p 3> I, all the components of the correlation tensors drop to zero, see 
FigGl This decrease is exponentially fast for the xz and zx components. For the 
other components, the asymptotic behaviour for large p simply scales like l/p 3 and 
does not depend on I anymore. This follows from the large argument asymptotics 
M 3 / 2 ~ 2^ an d -^5 /2 ~ tJ= ■ Note that in this case, we recover an algebraic 
decay similar to the local medium case given in Eq.0. 

To summarize this section, we have shown that when we take into account the 
non-local nature of matter by introducing a locality length I for the sources of the 
field, the correlation length is about I when the distance to the interface z < I. 
In this regime, the energy density saturates to a value given by the electrostatic 
energy density expression taken in z = I. 

4 Nonlocal plasma 

We consider in this section another simple situation where the field correlation 
function can be calculated fairly easily. Restricting ourselves to a field in thermal 
equilibrium between field and surface, we use directly the fluctuation-dissipation 
theorem for the field. The relevant information is thus encoded in the electric Green 
tensor (i.e., the field's response function). The Green tensor contains a part due to 
the reflection from the surface that is actually dominating in the sub-wavelength 
distance regime we are interested in. We first review the corresponding reflection 
coefficients for an electron plasma, taking into account the finite response time 
of the electrons and their scattering. These two effects make the plasma behave 
like a nonlocal medium and give rise to the so-called anomalous skin effect. We 
then discuss the large-wavevector asymptotics of the reflection coefficients and 
the corresponding limits on the spatial coherence function. It turns out that the 
scattering mean free path is one key quantity that limits the coherence length at 
short observation distances. 
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4.1 Reflection coefficients 

We focus here on the electronic contribution to the dielectric function and describe 
the background ions, interband absorption etc. by a local permittivity Eb to avoid 
unnecessary complications. As is well known 1 29 1, the dielectric function of a bulk 
plasma is actually a tensor with two distinct spatial Fourier coefficients, a "longi- 
tudinal" si (q) and a "transverse" e t (q) where q is the modulus of the wavevector. 
(As before, we suppress the frequency dependence for simplicity.) The fields out- 
side the metal surface are characterized by the reflection coefficients r StP (K) that 
depend only on the magnitude K = |K| of the incident wavevector projected onto 
the interface. Out of the two polarizations s and p, we need in the following only 
r p (K) in the (non-retarded) regime K 3> oj/c. This coefficient is given, e.g., in 
the review paper by Ford and Weber 1 30 1 : 

pV ' 1/Z P (K) + 1 

We use here a dimensionless surface impedance Z p (K) that reads in the non- 
retarded limit 

oo 

^ <? =* 2 +ft (12) 



it involves the longitudinal dielectric function only for which we take the Lindhard 
formula 1301311 

3f? 2 u 2 Mu) 
ei{q)=e b + — — --jn (13) 

LO + \V U! + WJl (U) 

V + W 1 ,1 AS 

u = = — (14) 

qvF qi 

The plasma frequency is given by fP = ne 2 /(meo) with n, ~e, m the electron 
density, charge, and mass, respectively. 

From the nonlocal permittivity dl3H15> . two characteristic length scales can be 
read off: the mean free path Z m f p = vp/v and Vf/u>, the maximum distance over 
which an electron at the Fermi energy can move ballistically during one period 
of the applied electric field. In the following, we use the complex length I = 
Vf/{u + if) defined in dl4i to simplify the notation. 

The Lindhard formula, Eqs. dl3H15> . is based on a semiclassical description 
of the electron gas (classical particles with Fermi statistics) with a damping rate 
v and a velocity vf at the Fermi energy. This description is valid as long as q is 
much smaller than the Fermi wave vector Hf = mvp /ft- Our model thus applies 
reasonably well to a "clean metal" where the mean free path is much longer than 
the Fermi wavelength, and to distances above 1/kp (typically a few A). Ref . 1 30 1 
gives a more general dielectric function that covers the regime q > kp as well. 
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4.2 Coherence function 

The fluctuation-dissipation theorem for the electric field, combined with the Green 
tensor describing the reflection from a planar surface, gives the following integral 
representation for the field's coherence function: 

oo 

E ij {v u v 2 ) = ^0{^T) / ^ £ Cf{Kp)^ T -^^- (16) 

with 7i = (u 2 /c 2 - K 2 ) 1 / 2 (Im7i > 0). For more details, see for example |32 
1331 . We have omitted the free-space part of the Green tensor that gives the same 
result as for the blackbody field. This part actually becomes negligible compared 
to the surface part given here if we focus on the sub-wavelength regime, z\ = z 2 = 
z«A: the integration domain ui/c < K < oo (which is absent in the free-space 
field) then makes the dominant contribution to the integral. 

The tensors (Kp) in d!6i depend on the lateral (signed) distance p = 
x\ — X2, as introduced after Eq.@. In p-polarization, it is given by 

7^2 2 / Jo — <h 2Ji\ 
C {p \K P ) = — r J + J 2 , (17) 
luJ \ -2Jx 2 J J 

involving the Bessel functions J„ = J n (Kp), n — 0, 1, 2. A similar expression 
applies in s-polarization. We can focus, for short vertical distances, on the range 
oj/c <C K, expand the reflection coefficients and find that \r s \ \r p \; hence, 
the s-polarization is neglected in the following. This also justifies our taking the 
non-retarded limit of the reflection coefficient (II Q . To the same accuracy, we ap- 
proximate 71 w i\K\. Finally, the correlation tensor becomes 

00 

£ 13 {v 1 ,v 2 )=^uj0{uj,T) I ^e- 2Kz £ C^iKp^mr^K). (18) 

M=s,P 

We anticipate from the integral representation (1181 that the wave-vector depen- 
dence of Im r p (K) determines the spatial coherence length: if K c is the scale on 
which lmr p (K) — > 0, we expect that the divergence of the energy density is 
smoothed out for z <C 1/ K c and that the lateral coherence length remains finite: 

4oh - l/K c for z < l/K c . 

4.3 Local medium 

Let us illustrate first how the Lindhard reflection coefficient reduces to its local 
form (the Fresnel formula). If the g-dependence of £i(q) can be neglected, writing 
£l(q) — > eioo the surface impedance H2\ integrates to Z p — > l/e\ oc . Eq.(^3 then 
recovers the reflection coefficient for electrostatic images, r p — > (ei oc — 1)/ (ei oc + 
1 ) which is the large K limit of the Fresnel formula for transverse magnetic (TM or 
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p) polarization. The integration of the Bessel functions and exponentials over K in 
can be carried out, and we get Eq.@ with its unphysical 1 /z 3 divergence. 
The same divergence would be obtained here from the background permittivity 
Eh that we assume local. To focus on the nonlocal contribution from the electron 
plasma, we consider the regime where E b is real so that the leading-order, local 
contribution analogous to Eq.(@J vanishes. 



4.4 Nonlocal reflection coefficient 



To get a qualitative insight into the impact of nonlocality, we perform an asymp- 
totic analysis of the dielectric function (I13H15> : 



£b - 
£b ( 1 



n 2 



lo{uj + iv) 

1 \ iC 



i + (gef 



11' 

3lJ 



q 2 A 2 



< 1 
> 1 



(19) 



where 



A = 



^/e b v 2 F /(3f2 2 



(20) 



is the Thomas-Fermi length that provides another length scale, and we use the 
notation C = 3Tru>f2 2 /v F . We recall that I is the complex characteristic length 
defined in dl4> . Note that for small q, we recover the usual, local Drude expression 
for an electron plasma 

n 2 

£loc = £& 7 -r^- (21) 

0J(U! + IV) 

At large q, one gets the dielectric function for Thomas-Fermi screening |29 1 with 
a screening length on the order of vf / & plus an imaginary correction. 

From the integral (I12> for the surface impedance, we find that the typical 
wavenumber is of the order of q > K. Hence the two limits quoted above translate 
into the following asymptotics of the reflection coefficient, after performing the 
integrations, 



Im r p (K) 



Im 



1 



£lo 



1' 



CKA 4 g{KA) 



3e 



b l + K/{e b ^K 2 + l/A 2 ) 



\Kl\ < 1, 
\Kl\ > 1. 



(22) 



The dimensionless function g(KA) is the integral 

oo 

g(KA) 1 (1/ 



^(KA) 2 + t 2 [(KA) 2 + 1 + t 2 f 



(23) 
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Fig. 3 Reflection coefficient Imr p (if) vs. the normalized wave vector Kvf/O. Dashed 
lines: asymptotical formulas <22> , Symbols: numerical calculation. The arrows mark, from 
left to right, the characteristic scales ui/c, l/\£\ and 1/A. Chosen parameters: fi/v = 192, 
c/vf = 148, (vf/(2 — 0.84 A), taken from the Drude description of aluminium 1291 . 
We take ej — 2 to model the contribution of bound electrons. Left panel: u) = 0.8 v or 
A = 19 fim. Right panel: cj = 0.55 fl (A = 140 nm), near the large- AT asymptote of the 
surface plasmon resonance in the local approximation (given by ei oc + 1 = 0). 



This can be evaluated in closed, but barely instructive form involving a hypergeo- 
metric function; its limiting behaviour is 

g(KA) a \n{l/KA) + In 2 - i for KA<. 1, 

sCKjI) = 1{KA)- 4 foxKA > 1. ( ' 

The first line applies to the intermediate case l/\£\ <C K ^ 1/yl, the second one to 
the regime A' > 1/|^|. In both cases, Eq.(|23 implies that \lmr p (K)\ < 1. 

The reflection coefficient is plotted in Fig[3] where the asymptotic expres- 
sions i22\ are represented as dashed lines. We find good agreement outside the 
crossover range K\£\ ~ 1. In the frequency range of the anomalous skin effect, 
ijj r~j v (left panel, A = 19 fim in the infrared), the nonlocal plasma shows an in- 
creased lmr p (K), with a cutoff occurring beyond K c ~ 1/yl [see Ea.( l24H . This 
effect is well known 1 30 1 and is related to the enhanced spontaneous emission rate 
for a nonlocal metallic surface that was recently pointed out |34|. The reflection 
loss remains small in absolute numbers because of the large conductivity of the 
material. The opposite behaviour is found near the (local, non-retarded) surface 
plasmon resonance (right panel, A = 140 nm in the far UV): Im r p (K) decreases 
from its local value, with a weakly resonant feature emerging around K ~ 1/|^|. 

From these plots, we observe that the characteristic wave vector scale K c 
strongly depends on the frequency range. An upper limit is set by 1/yl, involv- 
ing the Thomas-Fermi screening length, but significant changes already occur on 
the scale l/\£\. The characteristic distance below which non-local effects become 
manifest, is thus given by the largest of \£\ and A. This is typically \£\, since in 
order of magnitude, |^|/yl ~ Q/\oj + iv\ which is much larger than unity for 
good conductors up to the visible domain. At frequencies smaller (larger) than the 
damping rate v, the mean free path Z m f p (the "ballistic amplitude" Vp/iv): sets the 
scale for nonlocal effects, respectively. 
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Fig. 4 Electric energy density, normalized to its far-field value, vs. normalized distance 
fiz/vF- Dashed line: local dielectric. Solid line: numerical calculation (left: ui — 0.8 v\ 
right: lu — 0.55 J?; other parameters as in Fig|3J. The arrows mark, from left to right, the 
characteristic distances A, and A = 2-nc/uj. 



We note that for typical metals, the Thomas-Fermi scale A does not differ 
much from the Fermi wavelength l/fcj?. The asymptotics derived above within 
the semiclassical Lindhard model Jl 3I > is therefore only qualitatively valid at short 
distances (large wavevectors). 



4.5 Energy density and lateral coherence 



The numerical calculation of the correlation function Ea.( ll8t can be done effi- 
ciently using a numerical interpolation of Im r p (K) that we continue for large and 
small K using the asymptotics derived above. 

We plot in Fig|4]ihe electric energy density as a function of distance, for the 
same two frequencies as in FigO Deviations from the local approximation (dashed 
line) occur at distances smaller than enhancement at low frequencies (a; ~ v, 
left panel), suppression near the surface plasmon resonance (right panel), which is 
consistent with the trends seen in Fig|3] A saturation at small distances is also vis- 
ible, although it occurs for fairly small values of Qz/vp (where the semi-classical 
Lindhard function is in practice no longer valid). We note also that for z > A, the 
plots are only qualitative since the calculation does not take into account retarda- 
tion. 

Finally, we illustrate the finiteness of the coherence length as the distance of 
observation enters the nonlocal regime. We plot in Fig|5]the zz-component of the 
normalized coherence tensor (0, as a function of the lateral separation pjz. In 
the local regime, one gets a universal curve independent of the distance (dashed 
line). This is no longer true near a nonlocal metal: when Thomas-Fermi screening 
sets in (z < A), the coherence function departs from its local limit, its width (the 
coherence length) becoming much larger than z. 
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Fig. 5 Normalized degree of spatial coherence for z-polarized fields, probed at a lat- 
eral separation p. The numbers on the curves (solid lines) give the normalized distance 
Qz/vf = 100,10,1,0.1, with the normalized Thomas-Fermi screening length being 
QA/vf = (eb/3) 1 ^ 2 w 0.8. Dashed line: result for a local dielectric in the near-field 
limit z<A, taken from Eq.j4}. The chosen parameters are those of Fig[3] right panel. 

5 Concluding remarks 

We have discussed in this paper the impact of a nonlocal dielectric response on the 
spatial coherence of thermal electromagnetic near fields above a planar surface. 
Using two different models to describe the nonlocal response, we have shown that 
when the sources of the field have a finite correlation length, this length sets the 
minimum scale for the coherence length of the field as well. This behaviour is 
qualitatively similar to what we found previously when investigating the contribu- 
tion of thermally excited surface plasmons where coherence length and plasmon 
propagation length coincide 1 5 1 . We have thus provided semi-quantitative evidence 
for the impact of nonlocality that has been conjectured already by Rytov's group 



The calculation for an electron plasma model highlights, on the one hand, the 
crucial role played by Thomas-Fermi screening, that sets the minimum coherence 
length. On the other hand, significant deviations from the local description already 
occur at scales below the electron mean free path (Fig[3]and Fig|4j, although these 
are not accompanied by an increase in spatial coherence. 

Our calculations can be improved taking into account quantum effects in the 
Lindhard dielectric function |30|, which will lead to quantitative changes at short 
distance. Indeed, for typical metals, the Thomas-Fermi screening length vf/^ and 
the Fermi wavelength 1/kp are fairly close |29|. A comparison to other models 
of nonlocal dielectric functions would be interesting as well. On the experimental 
side, it would be interesting to compare the recent data on heat transfer between a 
scanning tip and a surface 1 13 1 with a microscopic calculation along the lines used 
here. We also mention that in the context of the Casimir force, nonlocal surface 
impedances have been studied. The nonlocal correction is particularly relevant 
at finite temperature and large distances and leads to a behaviour of the Casimir 
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force that is qualitatively similar, even without absorption, to the local, lossy Drude 
model, see for example Refs.[35 36 1. Finally, it remains to study the impact of 
another property of real metals, the smooth rather than abrupt transition of the 
electron density profile at the surface: this can be described by effective surface 
displacements that depend on both polarization and wave vector, thus adding to 
the nonlocal effects considered here 1571 . 

We thank Remi Carminati and Jean- Jacques Greffet for discussion and Illarion Dorofeyev 
and Francesco Intravaia for helpful comments. C.H. acknowledges support from the Euro- 
pean Commission (network FASTNet and projects ACQP and QUELE). 



A Appendix 

Les us consider the Green tensor relating an electric current in a local medium 2 
(z' < 0) to the electric field in medium 1 (z > 0) that we take as vacuum (si — 1). 
This tensor can be written as an expansion in plane waves (Weyl expansion) 

G,,(r,r') = '-J ^ gi^e^'-'^+^-^e^'e-^'' , (25) 

where K = (k x , k y ) is the wave vector component parallel to the interface. The 7$ 
are the z-components of the wave vector: 7? = e^w/c) 2 — K 2 . In the notation of 
Ref.Q, 

9ij( 

fl=S,p 

The polarization vectors for the s and p polarization are 



= e< 2) = K x e z (27) 

m Ki - 7K 

= t 1 — (28) 

P UJ/C 

z — 7 



where K is the unit vector parallel to K. The i 21 are the Fresnel transmission 
coefficients between media 2 and 1 : 

e = ^, tf = J2^L. ( 30) 

7i+72 1 £271 + 72 
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